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Abstract. A great deal of recent data on event-by-event fluctuation and correlation 
measurements has been released by several experiments at the SPS and RHIC. 
Recent results on charge fluctuations, balance functions in pseudorapidity, and average 
transverse momentum fluctuations will be reviewed. The results will be compared to 
various model predictions after examining contributions to each observable from known 
physics processes. 



1. Introduction 

Over the past several months, many results on event-by-event fluctuation and correlation 
observables have been made available by experiments at the SPS and at RHIC. This 
article will review the most recent results on three of these observables: net charge 
fluctuations, balance functions in pseudorapidity, and average transverse momentum 
fluctuations. Each of these observables is designed to probe a specific property of 
the system: net charge fluctuations are sensitive to the charge distribution within the 
collision volume, balance functions are sensitive to the hadronization time of the system, 
and average pt fluctuations are sensitive to critical fluctuations of the temperature of 
the system which may be present during a phase transition from hadronic matter to 
a Quark-Gluon Plasma (QGP). However, as with any measurement, the results can 
be contaminated by known processes such as charge conservation, resonance decays, 
elliptic flow, or even hard scattering processes. Therefore, the contribution of any known 
processes must be thoroughly examined prior to formulating a reliable interpretation of 
the results. Detailed background information on the theory behind various event-by- 
event observables can be found in other recent reviews pQ. 



2. Net Charge Fluctuations 

A promising signature for the presence of a QGP is the observation of event-by-event 
fluctuations in the net charge of produced particles. It is hypothesized that fluctuations 
in the net charge would be significantly reduced in a QGP scenario. This is due to 
the fact that in a QGP, the fractional electric charges of the quarks are more evenly 
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Table 1. The expected magnitude of various charge fluctuation measures within 
specific physical scenarios for a small sample of the total number of produced charged 
particles, ^ yn values are listed within both the STAR and PHENIX acceptances. 



Scenario v(Q) v + -, dyn (STAR) v+-, dyn (PHENIX) </> q 



Independent Emission 1.0 0.0 



0.0 



0.0 



Resonance Gas 0.75 -0.0013 

Quark Coalescence 0.83 -0.0008 
Quark-Gluon Plasma 0.25 -0.0038 



-0.006 
-0.004 
-0.019 



-0.125 
-0.084 
-0.375 



spread throughout the QGP volume than the unit electric charges of the hadrons in a 
hadronic gas volume. As a result, the fluctuations of the net charge in a given region of 
phase-space would be significantly reduced in a QGP 

There has been much debate on the choice of the most appropriate observable to 
quote for the magnitude of net charge fluctuations. A standard is yet to be agreed 
upon in the literature, so PHENIX j3] quotes the normalized net charge variance, v (Q), 
STAR [Sj quotes a measure for dynamical fluctuations, u + ^.^yn, and NA49 [H] quotes 
the variable 4> q . The latter will be quoted here. Fortunately, all of the variables listed 
above can be related to each other with some very reasonable assumptions [Zj. 

The baseline for the measurement of net charge fluctuations is that of independent 
particle emission (where v (Q) = 1.0, u + _^ dyn = 0.0, and (p q = 0.0). Predictions for the 
magnitudes of these three fluctuation variables for a variety of physical scenarios are 
listed in Table 1. Keep in mind that these values are valid only for small samplings of 
the total number of produced charged particles since charge conservation can contribute 
to any observed fluctuations. Net charge fluctuation measurements can be corrected 
for charge conservation based upon the fraction of the total number of particles in the 
sample as follows [Hj: 



Figure 1 compiles the results of net charge fluctuation measurements by NA49 in 
v /5^7 = 17 GeV Pb+Pb collisions [6J, and by PHENIX 4] and STAR [5, in ^sjfij = 
130 GeV Au+Au collisions, all quoted in terms of the (fi q variable. Not included are 
new results from the CERES experiment p]- Superimposed on each plot are curves 
representing the predictions for independent particle emission, quark coalescence 0, 
a resonance gas [3], and a QGP [3], after correcting for charge conservation within 
the acceptance of each detector. At SPS energies, NA49 observes fluctuations that 
are most consistent with the independent particle emission scenario. The PHENIX 
measurement is also consistent with independent particle emission. However, for 
a large pseudorapidity window, STAR reports fluctuations that are most consistent 
with the prediction of the quark coalescence model. Differences in the PHENIX and 
STAR measurements may be attributed to the reduced sensitivity of the PHENIX 
measurement, which may be improved in the analysis of Au+Au data at a/s/vtv — 200 
GeV. Note that all of the measurements differ dramatically from the expectations for 
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the pure QGP scenario, although it may be premature to rule out the possibility of the 
existence of a QGP based upon these results until more detailed models are consulted. 
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Figure 1. Charge fluctuation measurements as a function of acceptance for a) -^/s^vjv 
= 17 GcV Pb+Pb collisions from NA49 6 1, b) ^/JWn = 130 GeV Au+Au collisions 
from PHENIX gj, and c) ^/IWn = 130 GeV Au+Au collisions from STAR 5 . The 
curves represent expectations for independent particle emission, a resonance gas, quark 
coalescence, and a QGP after correction for charge conservation effects. 



3. Balance Functions 

The balance function analysis method is designed to determine whether hadronization 
from a system occurs early (< 1 fm/c) or at a later time ^0]- The method is based 
upon the fact that charge is conserved locally due to produced pairs that are initially 
correlated in coordinate space. Concentrating on balance functions in pseudorapidity, 
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the balance function is denned as: 



B(Ar)) 




(2) 



Here, (A?]) refers to a histogram of \rj(h + ) — rj(h~)\ for all possible pairs of the 

subscripted charge signs within a single event, summed over all events. N + and AL are 
the total number of positive and negative charged particles in the event, respectively. 
When examining the pseudorapidity of produced pairs if hadronization occurs at an 
early time, the pairs have more time to separate in 7] due to a combination of the 
expansion of the system and rescattering. Therefore, if there is a delayed hadronization 
scenario (as might be expected in central collisions), the analysis would yield a narrower 
balance function than the early hadronization scenario (as might be expected in more 
peripheral collisions) . 

Balance functions in pseudorapidity have been measured by STAR ^TJ and 
NA49 |12|. The behavior of the balance functions are studied as a function of centrality 
by plotting the widths of the balance functions in pseudorapidity (< At] >) obtained 
from a Gaussian fit or from a weighted mean in Arj. A reduction of the magnitude 
of B(Arj) for small Arj is observed due to the HBT effect, but it has little effect on 
the extracted widths. Figure 2 compares the balance function widths as a function of 
centrality at both SPS and RHIC energies. In both cases, there is a significant narrowing 
of the balance function observed for more central collisions. STAR has shown that the 
widths in the most central collisions can be reproduced by a thermal model, but only 
if an unrealistically high value of the transverse flow velocity, (3t = 0.77c is used with a 
low temperature of 105 MeV [TT]. Taking the interpretation at face value, this behavior 
is indicative of a late hadronization scenario. 

However, contributions from other known sources must first be considered when 
interpreting the results. An obvious contributor to the width of balance functions are 
neutral resonance decays, which will register in the balance function as produced particle 
pairs that occur at a late time, serving to broaden the function. The magnitude of 
resonance decay contributions using a thermal model baseline with (3t ~ 0.5c has been 
investigated by P. Bozek et al. [13]. The results of the calculation can reproduce widths 
as a function of centrality that are consistent with the STAR data. The behavior of 
balance functions as a function of centrality has also been investigated in the context 
of the quark coalescence model [H|. This model can also reproduce the centrality- 
dependence and magnitude of the widths of the STAR data with (3t ~ 0.5c, but without 
the inclusion of resonance decay contributions. 

4. Average px Fluctuations 

Event-by-event fluctuations in the average transverse momentum may be sensitive to 
temperature fluctuations near the QCD tri-critical point fH]- There have been many 
recent results from the SPS (NA49 [13 E2] and CERES [HI El) and RHIC (PHENIX 
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Figure 2. The widths of balance functions in pseudorapidity (< A77 >) as a function 
of centrality for a) ^snn = 17 GeV Pb+Pb collisions from NA49 and b) ^snn 
— 200 GeV Au+Au collisions from STAR JS]. In both cases, the balance functions 
are narrower for more central collisions. 



|19| I2U] and STAR |2UEI]) experiments. All experiments report a positive signal of non- 
random fluctuations, each with an apparent maximum in mid-peripheral collisions, as 
shown in Figure 3. As with the charge fluctuation measurements, a standard observable 
has not been agreed upon in the literature, and each experiment quotes a different 
variable to describe the magnitude of the fluctuations: <fi PT , F PT , cr^ ynamic , and S PT . 
Fortunately, the variables can all be related with very reasonable assumptions [22] • All 
of the variables measure the difference in the width of the event-by-event average pt 
distribution, M PT , from a random expectation, which can be expressed as a Gamma 
distribution that is directly calculable using the inclusive pr and number of particles 
per event, N, distributions [25] . 

Since a positive signal is observed and confirmed, contributions to the signal from 
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known sources should first be studied. The various analyses have addressed this 
with respect to contributions from HBT, resonance decays, and elliptic flow. The 
contributions from all of these sources have been found to be very small or negligible, 
each confirmed by independent analyses. Another possible source to the signal is that 
from hard scattering processes, which would be expected to produce correlations in 
with a stronger contribution at higher p?. 111 this context, both PHENIX and CERES 
observe that the majority of the contribution to the non-random fluctuations occur at 
high px- PHENIX observes a large 243% increase in the signal when the maximum of 
the px range over which the average p? is calculated increases from 1.0 to 2.0 GeV/c. 
This very large increase cannot be attributed to number fluctuations since the number 
of particles increases by less than 15% over that range. 

In order to investigate possible contributions due to jet production, PHENIX has 
applied a two-component model consisting of 1) a simulation of soft processes by 
reproducing the inclusive px and N distributions as a function of centrality, and 2) 
a simulation of hard scattering processes by embedding PYTHIA jet events at a given 
rate per produced particle in step 1 (while strictly conserving the N distribution so as 
not to introduce unnecessary number fluctuations). The simulation is performed as a 
function of centrality within the detector acceptance in two modes: a) with the PYTHIA 
event embedding rate kept constant for all centralities, and b) with the PYTHIA event 
embedding rate scaled by the nuclear modification factor, Raa, for each centrality bin. 
The result of the simulation superimposed with the PHENIX data is shown in the 
bottom plot of Figure 3. The simulation with the -R^-scaled rate agrees with the 
centrality-dependent (and also the pr-dependent jTH]) data trends remarkably well. 
Within this simulation, the decrease of fluctuations in the most peripheral collisions 
is attributed to the fact that the signal begins to compete with the magnitude of the 
number fluctuations, while most of the decrease of fluctuations in the most central 
collisions can be attributed to the onset of jet suppression. When the simulation is 
performed at ^/sjviv — 130 GeV within the STAR acceptance with the same PYTHIA 
event embedding rate that reproduces the PHENIX data, agreement with the STAR 
more central data is achieved within errors. The simulation predicts a 15-20% increase 
in the fluctuation signal from STAR at ^/snn = 200 GeV. Although the jet production 
cross section is much reduced at SPS energies, it is non-zero, and there may also be 
a residual jet contribution present. The identical simulation performed at ^/snn = 17 
GeV within the NA49 acceptance can also reproduce the trend of the NA49 data as a 
function of centrality. 

A promising method in which to directly compare the fluctuation measurements 
from all four experiments is to estimate the magnitude of any residual event-by- 
event temperature fluctuations present using the prescription outlined in [21]. The 
most central fluctuation values thus correspond to temperature fluctuations of 1.8% in 
PHENIX, 1.7% in STAR, 1.3% in CERES, and 0.6% in NA49. The residual temperature 
fluctuations are small and do not significantly increase between SPS and RHIC energies. 

Alternative interpretations of the observed average pt signal have been presented. 
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One such interpretation is that the signal may be due to the onset of thermal equilibrium 
within the system |22j . This model is successful in reproducing both the centrality- 
dependence of the fluctuation signal and the inclusive < pt > distributions from STAR 
and PHENIX. 




Figure 3. Average pr fluctuations as a function of centrality for a) ■ s /snn = 17 GeV 
Pb+Pb collisions from NA49 [T2], b) ^/JWn = 130 GeV Au+Au collisions from STAR 
IHl, and c) ^/IWn = 200 GeV Au+Au collisions from PHENIX [THj. The error bars 
include statistical and systematic errors. The curves are the results of a model that 
estimates the relative contribution due to jet production as a function of centrality 
[P3| for a fixed PYTHIA event embedding rate (dashed lines) and for an -R^-scaled 
embedding rate (solid lines). Uncertainties on the simulation are 0.3 MeV/c for NA49, 
2.1 McV/c for STAR, and 1.0 MeV/c for PHENIX. 
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5. Conclusion 



There has been nearly a doubling in the quantity of preliminary or published results 
on event-by-event fluctuation and correlation observables over the past year. Long 
thought of as the exclusive domain of large acceptance detectors, the scope of the 
results are demonstrating that event-by-event physics is accessible to smaller acceptance 
experiments with only a small sacrifice in sensitivity. The development of event-by-event 
analysis techniques is continuing to mature at a rapid pace, including extensions to 
multi-dimensional correlations and methods to relate fluctuation observables directly to 
correlation functions. It is becoming clear from what we have learned so far that these 
methods can provide additional tools with which to study the underlying dynamics of 
the collision system. 
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